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Abstract— Due to the existing variability in the calculation of the bearing 
capacity factors, specifically when obtaining the Ny value of the bearing 
capacity, which is estimated by a large number of equations that have 
of limit 
characteristics, empirical, semi-empirical, among others; thus existing a 


been obtained by approaches equilibrium, method of 
great uncertainty, therefore, this article took a total of 55 equations that 
compute this value to carry out a probabilistic study where Ny values were 
obtained for friction angles of 5°, 10°, 15°, 20°, 25°, 30° and 35°, being the 
variability obtained being directly proportional to friction angle. Also, it 
was found that the data have a good correlation with the Lognormal 
distribution, according to the Kolgomorov-Smirnov and Anderson-Darling 
tests, in this context, Lognormal cumulative probability distribution curves 
of the data under study were elaborated, thus allowing not to opt for a 
single or deterministic value but a probable or acceptable value 
according to the size of a project with a conservative range. 





I. INTRODUCTION 


The soil plays an important role as a support base for 
any construction, the gravitational load transmissions of 


Where q, is the ultimate bearing capacity of the soil, c is 
the cohesion, q is the overburden pressure, y is the unit 


weight of soil, B is the width of footing, N., N, and 


any structure go to this outcrop, where its load capacity 


indicates to the specialist a foundation geometry to choose 
from, and it is precisely that, in shallow foundations, 
bearing capacity generally governs the design process 
(Shahin & Cheung, 2011), representing a primary step for 
shallow foundation stability evaluations (Dewaikar et al., 


2008). 


For decades it has been known that bearing capacity 
has been studied and it has always been noted that it 
depends on certain factors; for example, without reaching 
various differentiated cases Terzaghi (1943) proposed an 
equation well known to date for the calculation of the 
bearing capacity considering the case of a strip footing, 


which is detailed below. 


1 
q, -cN, +qN, real 
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N, are the bearing capacity factors which depend on the 


friction angle ø of the soil. 


From equation 1, is particularly see that the values of 
N, and N, were made known by Prandtl in 1921 and 


Reissner in 1924; however, there is a great dissertation on 
the value N, since this is theoretically less precise than 


the other two terms (Griffiths, 1982), so there is a wide 
controversy considered over the theoretical values of N, 


(Mahmood, 2018) thus leading to great variability in the 
objective of achieving the value of the bearing capacity by 
various methods (Padmini, 2018; Shill & Hoque, 2015; 
Sieffert & Bay-Gress, 2000), is so the bearing capacity 
must be better understood using new parametric and 
(D) numerical analyzes (Sieffert & Bay-Gress, 2000), in this 
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sense, today there are different investigations based on 
analytical model approaches, semi-empirical models, 
empirical models, finite difference models, upper limit and 
lower limit models, finite element models, etc. that try to 
determine this geotechnical parameter (Motra et al., 2016). 
As described, there are several proposed formulas, but no 
is totally accurate (Nguyen et al., 2016), generating as 
consequences discrepancies in the results (Ty et al., 2019). 


The great geotechnical variability of the soil is well 
known, existing a certain proportion of uncertainty in the 
same parametric measurement, in this case, the N 7 value. 
Now, in this context, it is known that there are many forms 


of uncertainty as in JCGM (2008), 
geotechnical uncertainty is recognized to be primarily 


framed where 
epistemic, bayesian, and belief-based (Christian, 2004), of 
which this article will try to address the epistemic in a 
probabilistic way related to 55 formulas taken from Motra 
et al. (2016) that obtain the N, value, thus trying to 


improve the knowledge of this parameter, generating a 
reduction in epistemic uncertainties and therefore the total 
variability of the estimated soil design properties are 
reduced, as mentioned in Cao et al. (2017). 


Characterization of geological uncertainty remains 
challenging (Juang et al., 2019), and it is particularly 
knowing that the model to choose is still the main source 
of uncertainty (Motra et al., 2016), due to geotechnical 
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properties vary spatially (Han et al., 2020; Juang et al., 
2019; Popescu et al., 2005). However it is known that 
probabilistic methods provide a powerful tool for dealing 
with uncertainty in engineering projects (Christian, 2004), 
it is thus that what this article raises is to deal in a 
probabilistic way the value N, of bearing capacity under 


different friction angles ø , considering the results through 


a Cumulative Distribution Function (CDF), in this context, 
it will be allowed not only to opt for a single deterministic 
value but by a range of probabilistic values. 


I. METHODS 
2.1. Model for prediction bearing capacity factors. 


Although equation 1 is used in practice to calculate the 
bearing capacity of shallow foundations, however, to 
obtain the value N,, influencing the bearing capacity 


there are a large number of methods. Table 1 shows 55 
models to compute the N, value, which have been taken 


from the compilation of Motra et al. (2016). 


Table 1: Models for prediction N, value adapted from Motra et al. (2016) 






































N° Author Method Model 
1 Terzaghi (1943) Limit equilibrium N,- tan? (z + £) exp (x tan $) t 2o tan (1 34¢) 
adt P [ 2| T ó T ó 
2 Taylor (1948) Limit equilibrium N, =| tan T + 3 exp (x tan $) —1.0 | tan n + 
3 pe ee Method of characteristics | N, =| 1.413tan* ae $ exp(z tan) -1.794 | tan(1.27¢) 
(1953) L 4 2 
m 
4 | Biarez et al. (1961) Equilibrium limit N, = Ls un [Z4 f jexp(r tan é)-10| tang 
5 Feda (1961) Empirical N, = ooterp{ £) (for à < 35°, ¢ in degree) 
Semi-empirical based on (m @ 
N, =| tan”) —+=— Jexp(z tang )—1.0 |tan(1.4 
Wl 03) limit equilibrium , | | 4 £) pene) | (144) 
7 Hu (1964) Fitted model, equilibrium | y _/1 901tan?| Z +É lexp(tan)+0.27 |tan (1.285) 
limit 4 2 
T 69 "m 
8 Krizek (1965) Empirical N,- 40-4 (for ø < 35°, à in degree) 
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Table 1: Models for prediction N, value adapted from Motra et al. (2016) (continued) 
































N? Author Method Model 
9 Booker (1969) Method of characteristics | N, = 0.1045exp (9.69) 
10 Hang and Fitted model, Memad of N = tan? m " $ exp(z tang) ETT (1.339) 
Christensen (1969) characteristics i 4 2 
i T 
- = a Semi-empirical model | N, = 2.0 an (z 4. £) exp(z tan )-1 o| tang 
| eee Fitted model, limit | y=! 1.752 tan?/ Z+ Ê lexp(a tang) 0.186 | tan(1.32¢) 
Beik (1970) equilibrium í 4 2 
inch- m 
3 = Semi-empirical N, =1.5 a (z + £) exp(z tan 9) -1 ol tan ó 
14 Davis and Booker | Fitted model, equilibrium N. =| tan? m.p exp( tan $)+2.33 En (1.3164) 
(1971) limit , 4 2 
15 | Chummar(1972) | Fitted model, semi- | y _| 7.17 tan?/ 24% \exp(ztang)+65.5 [tan (0.274) 
empirical 4 2 
. a (zo 
16 Vesic (1973) Method of characteristics | N, = 2.0 an (z + d exp(z tan ) -1 o| tang 
17 Chen (1975a) Upper bound limit N, =2.0| tan? Z£ exp(z tan ø)+1.0 |tan tan TE 
analysis 4 2 4 5 





18 Chen (1975b) 


Fitted model, upper 
bound limit analysis 


N 


y 


f .45 tan? (z + d exp(z tan ø) + 0:754] tan (1.41¢) 



































(1993a) 


19 Salenon et al. Fitted model, limit N, — tan Zy (À exp(z tan 0) -10 tan (1.4054) 

(1976) equilibrium L 4 2 
20 Craig and Pariti Fitted Vra limit N, -|222tan? zb exp(a tan $)+0.222 tang 

(1978) equilibrium L 4 2 
21 Spangler and Approximation Dm N. 21.1 tan? a + 9 exp( tan $) —1.0 | tan (1.39) 

Handy (1982) Terzaghi's mechanism ' 4 2 
Ingra and Baecher nae : _ : 
22 (1983) Statistical analysis N, = exp(0.1739 —1.646) (¢ in degree) 
23 Sone ann Method of characteristics | N, =| tan? f Z+ p exp(z tan) 1.0 | tan (1.3419) 
Restaino (1984) d 4 2 
Hettler and - T 115] 
A Gudehus (1988) P iA TUM [57 (0n 9?) l 1-0 
: 0.757 
a E Limit equilibrium | N, = ew[ ing 5152896 -3452] 
m 

26 Bolonauo tan Method of characteristics | N, — [ws (z F £ exp (x tan $) =] ol tan (1 56) 
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Table 1: Models for prediction N, value adapted from Motra et al. (2016) (continued) 







































































N° Author Method Model 
2; | BoltonandLau | Fitted model, method of | y |) o74tan?{ “+ Jexp(atan) +3.736 [tan (1.3679) 
(1993b) characteristics 4 2 
ae | annona | en | ue Í 2tan? (z + d exp(z'tan d) 1 E tan (1.4174) 
solution 4 2 
29 Zadroga (1994) Empirical model N, = 0.657 exp (0.1419) (9 in degree) 
30 Manoharan and Fitted model finite N, = et Zr $ Jexp( tan ó) 43.464 an (1.2799) 
Dasgupta (1995) element 4 2 
Fitted model from K ,, tang( Ky 0.34 
31 Bowles (1996) E N= 9| An -1.0 |, K,, =exp| 1.708 +3.287¢—-—— 
values 2 \cos¢ i Ing 
32 Frydman and Burd Fitted model, Hoe N, zt a M (À exp (7 tan ó) sio e (I. 4 ó) 
(1997) difference analysis 4 2 
33 oco Upper P oo N, =exp(0.66+5.11 tan d) tan ø 
analysis 
34 Paolucci and Pecker Fitted mödel, upper N, = nèl Za ó exp(z tan $) 41.0 tan (1.71 $) 
(1997) bound limit analysis 4 2 
1.5 
Danish Standards STE _l (az ¢ 
35 Association (1998) Empirical fitting N,- 4 [o (z + $ Jexp(a tan $) -10] cos 7 
36 | Soubra (1999) Fitted model, upper y 11 374tan?( Z +2 lexp(atan 4) 0.162 [tan (1.3436) 
bound analysis 4 2 
2 m 
"T Coduto (2001) Approximation from | .374tan* (z + ^ exp (x tan $) +1 o| tan ó 
T hi’ l N,= 
Peau TUQUE , 1.0+0.4sin (4.06) 
Solution based on Davis 
Poulos et al. (2001 N, =0.1054exp(9.6 
98. | Poulos etal (2001 and Booker (1971) , POLA) 
i a 
39 | Uenoetal. (2001) | Fittedmodel, method of | y _| tani] 18 lexp(tang)—1.0 |tan(1.436) 
characteristics 4 2 
Fitted model one, upper mg: £) 
4 W t al. (2001 N, «1.2 tan | —+=— Jexp(zrtan ø )+ 4.6 |tan(1.436 
? npe AON bound limit analysis á | (5 2 2 9) ( 9) 
Fitted model two, upper [ Hg: £) 
41 |W t al. (2001b N, =| 1.234tan*| —+— Jexp(z tan ¢)+4.151 |tan (1.394 
ADEM bound limit analysis 7? (5 2 P( 9) ( 9) 
= 7 n 
42 | Zhuetal.(2001a) | Case 1, limit equilibrium | N, =| 2.0tan’ (s + d exp(zrtan)--1 o (tang) 
[ cm 
43 Zhu et al. (2001b) Case 2, limit equilibrium | N, = 2.0 tan* (z + £ exp (x tan $) +1 ol tan (1 .07¢) 
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Table 1: Models for prediction N, adapted from Motra et al. (2016) (continued) 


















































N° Author Method Model 
Dewakirand Fitted model, limit r n 
44 ilibrium based N, 2|1.626tan? | —+ = |exp(z tan $)-- 2.019 |tan (1.373 
Mohapatra (2003) equili nem ased on | n ? p( $) | ( $) 
Terzaghi's model = 
45 irae Gey | PE eset mend o | ay =| gagn f. g exp(z tang) + 0.508 | tan (1.3524) 
characteristics L 4 2 
46 | Kumar (2003b) Fitted model, upper | yw 11. 379tan?| 742 exp( tan 9) — 0.461 | tan (1.3374) 
bound analysis L 4 2 
47 | Hjiajet al. (2005) | Lowerandupperbound | y _ exp| Z (1+37 tan 9) (tand)* 
analysis 6 
ak |. amia | PPE d method OE | a | dan” (z + £ jer tanø)-1.0 [tan (1.3384) 
characteristics 4 2 
m 
49 Smith (2005) Method of characteristics | N, =1.75 [ws (z + £) exp [ (0.752 + $) tan g] -1 o| tan ø 
so | Far and Kanzer, | Lowerandupperbound | y 11.012 tan?/ £ L9 laplenn a 
(2007) limit analysis L 4 2 
51 Lyamin et al. Lower and upper bound N, = etn’ E t $ exp (x tan $) — 0.6 |tan (1 33) 
(2007) Am 4 2 
! Fitted model, lower r ] 
r r mT 
52 Kunar and kaiti bound finite element N, =| tan (4 4 £) exp(z tan ¢)—1.0 |tan(1.26¢) 
(2008) ] ] 4 2 
linear programming = 5 
Approximation model = = 
53 Salgado (2008) from N, values of N, =| tan £f Jexp( tan 4)-1.0 |tan (1.32¢) 
Martin (2005) 7 T 
54 Tanp and Yang Fatted model, upper N, =| tan? 2i ? exp(z tan 9) 1.0 | tan (1.3969) 
(2008) bound limit analysis L 4 2 | 
55 Jahanandish et al. pu) model, zero N. =| tan? m x (À exp (tan ó) 41.0 Itan (1.5) 
(2010) extension lines method L 4 2 | 




















As can be seen from the previous table, there is a great 
variety of models to obtain the N, value. This is due to 


the methodology used by the authors, the considerations 
and assumptions made such as the one usually used in 
which it is considered that the soil is a rigid material and 
that it generally fails due to shear (Meyerhof, 1951), and 
also due to the premise that soils may be homogeneous in 
terms of composition, they may not be homogeneous in 
terms of mechanical behavior (Uzielli et al., 2007). 


2.2. Variability of prediction models. 


It is known that geotechnics is not an exact science 


(Lacasse & Nadim, 1998) and that the specific 
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geotechnical parameters remain an open question in its 
analysis (Cao et al., 2016), this is clearly reflected in the 
large number of formulas presented in Table 1, that is why 
these formulas have been grouped for friction angles of 5°, 
10°, 15°, 20°, 25°, 30° y 35°; the value of up to 35° has been 
chosen because some formulas establish an approach of 
9 «35? e.g. Krizek (1965). The result of the variability of 


the N, values can be seen in the following figures. 
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Fig. 1: Availability of the N, values for each author 






5 Rm 


Q 
x 
tor © 


mq 
"1 
E 
P 
& 





v — —T 


E 
10 20 30 40 50 
Models for Ny 





Fig. 2: Variability of the N, values for the 55 models in 
the study 


From Fig. 2 shown, it is easy to check what is 
established in Padmini (2018) where the bearing capacity 
factor N, varies sharply with friction angle, especially 


when it is on the rise, thus Table 2 denotes both the 
maximum and minimum N, values for the studied 


database, where there is clearly a differentiation of values 
of up to 65.33 (specifically in ø 235?). Also, it is shown in 


Fig. 3 that the average values and the average values tlo 
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(standard deviation) for different friction angles, where 
most of the estimated models (ordered according to Table 
1) are within this range, reflecting the existence of 
similarity of values between model and model, but when 
wanting to consider a unique estimated value of N, , this 


becomes somewhat complicated. Although different 
countries already have predefined different methods used 
to estimate the N, factor as shown in Sieffert & Bay- 


Gress (2000), it may be that these values used in a 
deterministic way do not exactly define or characterize the 
soil profile for different geotechnical zones, as it is widely 
known that the variability of soil properties encountered in 
any project is related to the particular site and specific 
regional geology (Baecher & Christian, 2003), so it is 
convenient to see this problem through a statistical 
approach. Therefore, if we expand the methodology in a 
method of probabilities based on uncertainties, we will 
have the advantage of providing more complete and 
realistic information regarding the level of safety of design 
(Uzielli et al., 2007), finding not a true value, but an 
acceptable or conforming value, to the reality of the safety 
of any project. 


Table 2: Maximum and minimum values of N, 















































! Maximum and minimum values of N, 
$ ©) Maximum | Author | Minimum | Author 
5 1.81 pco a dem 
10 3.92 pco n Des 
15 6.62 peus US d 
20 11.02 $ p oar [48 den 
25 2231 pros is "en 
30 45.15 i ^ d 11.27 pe iy 
35 91.37 7 poy 26.04 pen 
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Fig. 3: Range of values *1lo for friction angles: a) 5°, b) 10°, c) 15°, d) 20°, e) 25°, f) 30° and g) 35° 


2.3. Probabilistic approach. 


The use of probability models and statistical methods 
to analyze data has become common practice in virtually 
all scientific disciplines (Devore, 2008), many studies that 
are based on probabilistic terms about geotechnics and 
soils have been developed in recent years (Al-Bittar & 
Soubra, 2013; Chen et al., 2012; Griffiths et al., 2002; 
Jung et al., 2008; Shahin & Cheung, 2011). If we see, for 
example, in a general approach, for engineering, the most 
common continuous probability distribution functions are 
the Exponential, Gamma, Beta, Uniform, Weibull, 
Rayleigh, Normal and Lognormal. On the other hand, in a 
slightly more particular approach, in geotechnics, the 
distributions used in the literature to model soil properties 
are the Lognormal, Gamma and Beta (Chen et al., 2012; 
Popescu et al., 2005). In this context, some examples are 
cited such as Vessia et al. (2009), who assumed in his 
study a Lognormal distribution, and Han et al. (2020) used 
the Beta, Gamma and Lognormal distributions. 


2.4. Probabilistic function. 


The probability function used in this article is the 
Lognormal. This distribution showed a good correlation 
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with the estimated N, values for the seven friction angles 


under study, which could be verified using both the 
Kolgomorov-Smirnov test (Massey, 1951) and Anderson- 
Darling test (Rahman et al., 2006) to ensure the goodness 
of fit of the distribution. The results are shown in Table 3. 


The Probability Density Function (PDF) of the 
Lognormal distribution is given as: 


Em "E :-) Q) 





1 
f(x; 4,0) = 
V270,,x 


On 
Also, the probability value associated with the 
Cumulative Distribution Function (CDF) is as follows: 
1 = 
P[X <x] -o| =H) GB) 
On 


Where 44, and o, are the mean and standard deviation of 


In(x) respectively. 


Equation 2 was used to match the probability 


distribution presented by the N, value at different friction 


angles. It was found that the variability of the data has a 
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good correlation with the Lognormal distribution as shown friction angles was plotted on the X-axis and the value on 
in Fig. 4, where the data of the N, values for the different the Y-axis Lognormal probability. 


Table 3: Tests for the estimation of N, probability Lognormal distribution 















































Testing 
itical Kol - 

¢ ; ee o naicdpdi Anderson-Darling test Em 

i significance Smirnov test Decision 

© level i Jn 

Critical T Critical "" 
Statistics Statistics 
value value 

5 0.05 0.17981 0.12815 2.5018 1.24130 The hypothesis should not be rejected 
10 0.05 0.17981 0.09717 2.5018 0.72867 The hypothesis should not be rejected 
15 0.05 0.17981 0.12073 2.5018 0.89589 The hypothesis should not be rejected 
20 0.05 0.17981 0.11905 2.5018 0.80133 The hypothesis should not be rejected 
25 0.05 0.17981 0.11198 2.5018 0.57770 The hypothesis should not be rejected 
30 0.05 0.17981 0.15610 2.5018 0.93521 The hypothesis should not be rejected 
35 0.05 0.17981 0.11978 2.5018 1.06750 The hypothesis should not be rejected 
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Fig. 4: Lognormal probability distribution curves for N, values for friction angles: a) 5°, b) 10°, c) 15%, d) 20°, e) 25°, f) 30° 
and g) 35° 
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III. RESULTS 


Based on the good correlation that exists between the 
Lognormal probability distribution, the data under study 
shown in the previous section and using equation 3 for the 
CDF, seven probability distribution curves were obtained 
for the friction angles 5°, 10°, 15°, 20° 25°, 30? y 35°, 
which are shown in Fig. 5. 














o 





Probability 
o o 











Fig. 5: Lognormal CDF of N, values for different friction 


angles 


IV. DISCUSSION 


The curves shown in Fig.5 are commonly used to 
express probabilities of damage under a certain level of 
demand (Ahmed & Soubra, 2012; Popescu et al., 2005) but 
they are also used to express probability levels of different 
parameters as in Vessia et al. (2009) who used it in the 
dataset of bearing capacity of strip foundation or Chen et 
al. (2012) who used it in undrained shear strength. In this 
context, these curves allow access not only to select a 
single value, that is, to a characteristic value of the 
parameter N,; instead, select a value based on 
probabilities, under certain reliability that the user 
considers relevant, thus allowing to reduce the uncertainty 
when choosing which appropriate equation best estimates 
this parameter. 


V. CONCLUSION 
In this article, the probability distribution of 55 
approaches by various authors was investigated to obtain 
the N, value of the bearing capacity under the friction 
angles: 5°, 10°, 15°, 20°, 25°, 30? y 35°. The results are 
shown as follows. 


There is great variability in obtaining the N, values 


through the 55 equations proposed by various authors, 
finding that this variation is proportional to the friction 
angle; that is, the greater the friction angle, the greater the 
variability of the estimated value of N, and vice versa. 
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Also, the N, values calculated for the friction angles 


mentioned above follow a suitable Lognormal probability 
distribution. 


Finally, faced with the problems encountered in 
engineering, specifically concerning soils, it is very 
common to use a deterministic practical approach; 
however, in parameters with great dissertation such as the 
N, value, key in the bearing capacity, it is not advisable 


to choose a unique value knowing that there are currently 
different methodologies to calculate, in that sense, by 
incorporating previously verified CDF, in obtaining this 
parameter, it was found that the results are shown using a 
probabilistic approach reduce bias and uncertainty to 
characterize the value of bearing capacity, generating in 
decision making opt not only for a deterministic value 
based on a single formula but to choose one or more 
probable or appropriate values to the conservative quality 
required in a project. 
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